Flexoelectric surface switching of bistable nematic devices 
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We report on a novel method of dynamically controlling the boundary conditions at the surface 
of a nematic liquid crystal using a surface flexoelectric effect. By moving the surface directors we 
show that one is able to manipulate defects which lie near the surface. This can be used to produce 
switching of a nematic liquid crystal device between two states with very similar free energies. This 
results in a bistable device that can retain either state with no applied voltage. Switching between 
the states occurs when the movement of the surface directors rotates those in the bulk which are 
then able to create or annihilate defects which lie near the surface of the device. 
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Defect structures in liquid crystals naturally evoke in- 
terest due to the appearance of complex textures easily 
visible by the naked eye. Early studies focused on clas- 
sifying the static properties of the defects and their in- 
teractions More recently the focus has moved to 
examining the dynamics of topological defects. 

In most practical applications of liquid crystals, such 
as traditional display devices, defects destroy the optical 
properties and are undesirable. However novel display 
designs, such as bistable displays or multidomain nemat- 
ics, exploit defect properties. Most attempts to control 
the defect motion have made use of bulk electric fields . 
In this Letter we will show how the flexoelectric effect can 
be used to control the surface alignment of a liquid crys- 
tal, thereby manipulating defects near the surface and 
switching the state of a liquid crystal display device. 

Liquid crystals are typically comprised of highly 
anisotropic, rod-shaped molecules. The nematic phase, 
which we shall concentrate on here, occurs when the 
molecules align parallel giving rise to long-range orien- 
tational order with the direction of alignment indicated 
by the so-called director field In a typical display 
device, the liquid crystal is confined between two plates 
a few microns apart. The director configuration on the 
plates is fixed. When an electric field is switched on the 
molecules align in the direction preferred by the field. Af- 
ter switching off the field, long-range elastic interactions 
ensure that the molecules reorient themselves in the di- 
rection preferred by the surfaces. These devices can be 
used as displays because different liquid crystal orienta- 
tions have different optical properties. The switching of 
such traditional liquid crystal devices is well understood. 

However there is now considerable interest in develop- 
ing bistable devices which can retain a memory of two 
distinct states, with different liquid crystal orientations, 
and hence different optical properties, even when the field 
is switched off [^|J^. It has been demonstrated, both by 
minimizing a Landau-de-Gennes free energy functional 
and by experiments [^j^, that it is possible to produce 
two (meta)stable states with different orientations of the 



director in the bulk by having different configurations of 
defects near the surface. However, the process of switch- 
ing between these two states is not understood. As the 
switching involves changes in the configuration of the de- 
fects near the surface, controlling this process requires 
understanding the dynamics of defects and how they can 
be manipulated. 

In this Letter we describe the physics behind the 
switching dynamics of a simple bistable nematic device. 
Our main conclusions are that the driving force which 
causes the switching is a surface flexoelectric effect. This 
rotates the surface directors which in turn pull round the 
bulk director field allowing the defects which lie near the 
surface of the device in one of the bistable configurations 
to be formed or annihilated. Once this is done switch- 
ing is essentially complete and the driving field can be 
removed. 

The main barrier to understanding the dynamics of liq- 
uid crystals is the complexity of the liquid crystal equa- 
tions of motion. The hydrodynamics of liquid crystals 
are usually described by the Ericksen-Leslie-Parodi equa- 
tions [||. However these are not sufficient to describe 
switching in bistable devices because they are restricted 
to an order parameter of constant magnitude. Defect mo- 
tion, which plays a major role in the bistable switching, 
is not included. 

Therefore it is necessary to consider a more general for- 
malism of the hydrodynamics in terms of a tensor order 
parameter, Q 0. Q evolves according to a convection- 
diffusion equation 



{dt + u- V)Q - S(W, Q) = FH 



(1) 



where u is the bulk fiuid velocity and F is a collective ro- 
tational diffusion constant. The term on the right-hand 
side of Eqn. describes the relaxation of the order pa- 
rameter towards the minimum of the free energy T , 



(2) 
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The order parameter distribution can be both rotated 
and stretched by flow gradients. This is described by the 
term on the left-hand side 

S(W, Q) = (CD + n){Q + 1/3) + (Q + I/3)(eD - O) 
-2e(Q + I/3)Tr(QW) (3) 

where D = (W + W^) /2 and O = (W - W^) /2 are the 
symmetric part and the anti-symmetric part respectively 
of the velocity gradient tensor Waf3 = dfjUa and ^ is re- 
lated to the aspect ratio of the molecules. The velocity 
field u obeys the continuity equation and a Navier-Stokes 
equation with a stress tensor generalized to describe the 
flow of liquid crystals. The details of the equations of 
motion can be found in reference |@]. We stress that the 
equations include backflow, allow variations in the magni- 
tude of the order parameter and follow the hydrodynam- 
ics of topological defects. To solve them we use a recent 
lattice Boltzmann approach which is proving particularly 
robust for problems which relate to complex fluids 

We consider a liquid crystal described by the Landau- 
de-Gennes free energy B 



T= [ dV{Fi,-FE+ Fd} + [ dS {Fa} , 

Jv JdV 



(4) 



where 



A 



Fe — -r^PaQapEp + {dpQal3){£fiEa + ef^Qa-fE^), 
izTT 

Fd = ^id^Qpxf, Fa^'^{Q^p^Q%f. (5) 



(Greek subscripts represent Cartesian directions and the 
usual summation over repeated indices is assumed.) The 
bulk free energy terms Fb describe a liquid crystal with 
a first-order, isotropic-nematic transition at 7 = 2.7 0. 
Contributions to the free energy which arise from an im- 
posed electric field E are included in Fe ■ Note that there 
are bulk terms which depend on Q and flexoelectric terms 
which couple to the derivatives of Q ||l^,|ll| . It is the lat- 
ter that we shall focus on in this paper. They arise in 
liquid crystals where molecules are asymmetric in shape 
and carry a dipole moment. Alignment of the dipoles by 
a field can then lead to a varying director field being en- 
ergetically favorable. Fd describes the elastic free energy 
within a one elastic constant approximation |^ and Fa is 
a surface free energy which fixes a preferred orientation 
Q° for the surface director field ||l2| . 

We work in two dimensions and consider the geome- 
try shown in Fig.^. The device width is and periodic 
boundary conditions are imposed in the y-direction. The 
surface pinning potential at a; = is periodic with the 
director preferring to lie in the (x, y) plane at an angle 
6*0 = (787r/180) sin(27ry/Lj,) to the x-axis. At x = Lx we 
take 6*0 = Ttt/ISO. 
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FIG. 1. The geometry used to demonstrate flexoelectric 
switching in a bistable nematic device. In particular note the 
orientation of the director on the surfaces of the device. The 
bold arrows show the direction in which the surface directors 
rotate when a field is applied in the negative a;-direction. 

Our aim is to show that in a device of this geometry 
switching can be driven by a surface flexoelectric effect 
and to demonstrate the kinetic pathway along which the 
switching proceeds. Therefore we put and 6/2 in the 
free energy expression equal to zero ||l3|| . Integrating 
by parts shows that the remaining fiexoelectric term can 
be rewritten as a surface contribution 



av 



dS {apEaQap} : 



(6) 



where a is the surface normal. In the geometry of Fig. |l| 
we will apply an electric field along x. Symmetry sug- 
gests, and we find in our simulations, that the director 
remains in the (x, y) plane. In this case, assuming uni- 
axiality (so that Qxx — 9^ (cos^ 9 — \ /3) where q is the 
related to the largest eigenvalue of Q), and ignoring elas- 
tic and bulk contributions the surface free energy density 
at X = is 



-60)^ +qhf,Ex{cosH-l/3) (7) 



where 9 is the angle the director makes with the x-axis. 
Minimising with respect to 9 gives, for small deviations. 



0-9a = ef,Ex/as sin 261. 



(8) 



Switching will be driven by a series of voltage pulses 

= 0, < t < io, (9) 

V^+Vo, to<t<to+ti, (10) 

V = 0, to + h <t <2ta+ti, (11) 

V^-Vo, 2to + ti <t <2ta + 2ti. (12) 

The first voltage pulse corresponds to a field in the nega- 
tive X direction. Then for < 9 < tt/2, 9 — 9o is negative; 
for —tt/2 < 9 < 0, 9 — 9o is positive. Therefore all the 
surface directors move towards the vertical as shown by 
the arrows in Fig.^. This will tend to favour vertical 
alignment, with the director field parallel to the x-axis, 
in the bulk. Conversely a field in the positive x-direction 
causes the directors to move away from 9 = favouring 
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a bulk state with the director field aligned diagonally as 
shown in Fig.||a. Thus the surface flexoelectric term can 
provide a driving force for the switching. 
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FIG. 2. Snapshots of the director configuration during 
switching: (a) no field, (b)-(e) switching from diagonal to 
vertical, (f) no field, (g)-(j) switching from vertical to diago- 
nal. Times are measured in simulation units. 

Simulations of the full Eqns. (0)-(||) displayed in Fig.|| 
show the path by which switching proceeds. At each time 
we display that portion of the device which lies near the 
X = surface. The voltage pulses cause the system to 
cycle between the states shown in Figs. ||a and |]f which 
are stable in the zero field portions of the voltage cy- 
cle ((^ and (^Tj) respectively). These have very similar 
free energies because of the imposed geometry. We find 
a free energy difference, measured as a fraction of the 
elastic energy, of 0.2%. The other frames in Fig. 2 show 
how the system moves between the bistable states. To 
switch from diagonal to vertical (Figs. 2a-2e) the surface 
deformation pulls the bulk directors towards the vertical. 
This provides the driving force for the pairs of defects at 
the surface to move together and annihilate. The switch 



from vertical to diagonal (Figs. 2f-2j) occurs in a similar 
way with the bulk rotating until there is sufficient driving 
force that the defects can be re-formed. 

Further insight into the process can be obtained by 
considering the free energy. This is shown in Fig. 3 where 
the bulk, elastic and total (elastic -I- bulk -I- surface -I- 
flexoelectric) free energies are plotted as a function of 
time. Upon switching the field there is a small spike in 
the free energy which quickly relaxes as the surface direc- 
tors move to their new positions. The switching itself is 
a far slower process - the elastic energy decreases slowly 
at the expense of the bulk free energy as the directors 
are pulled around by the surface until it is possible to 
push the defects out of or pull them into the device. De- 
fect removal or creation is marked by a maximum in the 
elastic free energy, which becomes more pronounced as 
the switching time becomes longer, together with a faster 
decrease in the bulk free energy. Once the elastic free en- 
ergy has passed this maximum switching will occur if the 
field is removed. 

The parameters used to obtain the results in Fig. 3 
were, in simulation units, A — 1.0, F = 0.5, 7 = 3.0, 
K = 2.0, i = 0.52, to = 2000, ti = 2000, L.^, = 63, Ly = 
24. The dimensionless ratio r = ejiE/as — 0.21 and 
the voltage was ramped up or down over 40 time steps. 
The simulation parameters can be mapped to physical 
values for comparison to prototype devices. However the 
mapping must be treated with some caution because the 
simulations are in two dimensions. We choose length, 
time and pressure scales so that — 1/im, pressure^ 
lO^Nm"^ viscosities ~ 0.04 - 0.1 Pa.s, k = 5 x 10""N 
and switching times ~ 10~'*s. 

r is a measure of the efficacy of the fiexoelectric effect 
in rotating the surface directors. For a given device this 
ratio is controlled by the applied voltage. As expected 
the switching time increases rapidly with decreasing r. 
As r decreases further switching fails. For example, for 
r — 0.05 and — 63, switching had not occurred after 
10ms. Moreover there was no discernible variation of the 
free energy with time suggesting that switching would 
not occur at a later time. 

It is also interesting to consider the effect of increasing 
Lx on switching times. The vertical to diagonal switch- 
ing time is little changed by doubling the device width, 
whereas the diagonal to vertical switching time increases 
substantially. This is because an increase in L^, with- 
out changing other parameters, favours the diagonal state 
over the vertical one: the difference in their free energies 
is now 5.8% of the elastic free energy. 

As K is decreased the switching time increases. This 
happens because the diffusion constant for the relaxation 
of the orientation of the director is proportional to k. For 
small K, switching fails. This occurs when the elastic cou- 
pling between the surface and bulk directors becomes too 
small for the bulk to be rotated by a surface perturbation. 



3 




5000 , 10000 



FIG. 3. Variation of the total, elastic and bulk free ener- 
gies (measured relative to the total free energy with no applied 
field) with time t. A field along —x is switched on at to = 2000 
and off at to + ti — 4000. A field along +x is switched on at 
2to +ti= 6000 and off at 2to + 2ti = 8000. Times, measured 
in simulation units, are comparable in Figs. 2 and 3. 

In the simulations it is possible to "turn off" back-flow 
effects, that is, to ensure that any movement of the di- 
rector field docs not induce a flow. This has little effect 
upon the switching time. Increasing the temperature to- 
wards that of the nematic-isotropic transition leads to 
a decrease in switching time because the defects have a 
lower pinning energy. 

The simple device described here has allowed us to ob- 
tain a fuller understanding of a possible switching mech- 
anism in bistable displays. Experiments on prototype 
devices have reported switching times somewhat longer 
than those calculated here. One likely reason for this is 
that we have tuned the free energy of the two metastable 
states to be very nearly equal. As we saw when changing 
the device width, the switching time is strongly depen- 
dent on this free energy difference. Moreover, in a real 
liquid crystal there will in general be a bulk field term 
in the free energy (e^ =/= 0). This may cause the switch- 
ing in one direction to be due to a Frederiks' transition 
and inhibit switching in the other direction. A change 
in boundary conditions could be used to compensate for 
this imbalance. Other differences in prototype devices 
are that the modulated director structure at the surface 
is created using normal alignment on a blazed surface 
rather than modulated alignment on a flat surface. Work 
to incorporate this additional physics is currently under- 



way. 

To conclude, we have shown that it is possible to model 
the switching mechanism in a bistable nematic device us- 
ing lattice Boltzmann simulations of the liquid crystal 
equations of motion. Thus we have been able to show 
that switching can be driven by a surface flexoelectric 
term in the free energy. Bistability occurs because the 
diagonal and vertical director alignments have very sim- 
ilar free energies. They remain stable because there is 
a free energy barrier between them. To promote switch- 
ing, defects need to be introduced to or expelled from 
the device: this is most easily accomplished if the defect 
pinning energy is small. 
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